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INTRODUCTION 

The  fourth  meeting  of  the  Mathematical  Model  Working  Group  was  held 
at  the  International  Center,  Stanford  Research  Institute  on  January  16,  17, 
and  18,  1973.  The  purpose  of  the  meeting  was  to  provide  the  opportunity 
for  detailed  Information  exchange  between  MMWG  members  on  matters  relating 
to  various  math  model  techniques  and  results. 

In  his  opening  remarks  Dr.  Sherman  (Working  Group  Chairman)  reviewed 
the  results  of  the  Transducer  and  Hydromechanics  Committee  program  develop- 
ment meeting  of  January  3,  4,  and  5.  During  the  THC  meeting  there  was 
considerable  discussion  and  Interest  in  the  need  for  noise  models  - 
particularly  as  applied  to  submarine  hull  structures. 

Present  at  the  SRI  meeting  were  the  following: 


C.  Sherman  (Chairman) 

- NUSC/NL 

R.  Radllnskl 

- NUSC/NL 

G.  Martin 

- NUC 

E.  Spurlock* 

- SRI 

R.  Smith 

- NUC 

J.  Young* 

- SRI 

J . Hunt 

- NUC 

R.  Kolesar 

- NSSC 

G.  Benthien 

- NUC 

C.  Campbell 

- ADL 

M.  Knit tel 

- NUC 

D.  Davison 

- NUC 

S.  Hanlsh 

- NRL 

P.  Rogers 

- NRL 

R.  Baler 

- NRL 

L.  Van  Buren 

- NRL 

R.  Dunham 

- NUSC/NL 

D.  Porter 

- NUSC/NL 

♦Part-time 


1 


MODEL  FOR  BAFFLING  OF  FLEXURAL  WAVES 


C.H.  Sherman 

Naval  Underwater  Systems  Center 


SUMMARY 


A very  simple  version  of  a model  for  evaluating  the  effectiveness  of 
baffle  structures  was  presented.  The  model  consists  of  an  infinite  plane 
elastic  plate  representing  the  hull  of  a ship  on  which  are  two  other 


plane  layers  representing  the  baffle.  The  layer  adjacent  to  the  hull 


layer  has  fluid  properties  while  the  other  layer  of  the  baffle  is  another 


elastic  plate.  Water  lies  beyond  this  second  plate.  The  fluid  layer  can 
be  given  a dispersive  velocity  and  other  frequency  dependent  properties 
to  make  it  represent  more  complicated  structures  such  as  a compliant  tube 
baffle.  The  outer  elastic  layer  is  probably  a necessary  part  of  the  baffle 
in  order  to  avoid  too  much  loss  of  hydrophone  sensitivity  when  hydrophones 
are  mounted  Just  outside  the  baffle. 

The  solutions  to  two  separate  problems  using  this  model  were  described 
The  first  is  the  calculation  of  the  effectiveness  of  the  baffle  layer  in 
reducing  the  noise  pressure  produced  by  flexural  waves  in  the  hull. 
Specifically,  the  ratio  of  the  pressure  at  the  outer  side  of  the  baffle 
to  the  pressure  on  the  hull  without  the  baffle  was  calculated  for  a 


specified  flexural  wave  in  the  hull 


The  second  problem  is  the  calculation  of  the  tendency  of  the  baffle 
layer  to  degrade  the  signal.  The  ratio  of  the  total  signal  pressure  on 
the  outside  of  the  baffle  to  the  total  signal  pressure  on  the  hull  without 
the  baffle  was  calculated  for  a specified  incident  sound  wave. 


■■^■■■■PPWRMMMPBMI 


The  ratio  of  the  two  ratios  described  above  is  the  signal-to-noise 
improvement  produced  by  the  baffle.  Although  this  model  is  oversimplified 
in  many  respects,  it  does  serve  as  an  example  of  the  type  of  calculation 
which  is  required  for  evaluating  different  baffle  concepts  and  other  aspects 
of  large  aperture,  hull-mounted  arrays. 
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FINITE  ELEMENT  ANALYSIS  OF  TOWED  ARRAY  MODULES 


R.R.  Smith 

Naval  Undersea  Center 


The  finite  element  method  ia  a numerical  technique  useful  in  the 
vibrational  analysis  of  complex  structures.  It  can  be  applied  to  towed 
arrays  giving  the  vibration  levels  (Including  hydrophone  output)  for  a 
variety  of  forcing  functions  (cable  strum,  turbulent  boundary  layer,  drogue 


whipping).  Ultimately,  the  finite  element  method  can  be  used  as  a design 


tool  to  predict  the  effect  of  various  construction  techniques  on  array 


self-noise  levels 


This  paper  gives  some  preliminary  results  using  the  finite  element 


method.  The  mathematical  model  is  shown  in  Figure  1.  It  consists  of  a 


module  33"  long  with  1.25"  0«D.  The  module  has  endcaps  on  both  ends  and 


two  hydrophones  freely  suspended  in  the  fluid.  The  left  end  is  vibrated 
axially  (to  simulate  the  axial  component  of  cable  strum)  and  the  right  end 


is  left  free.  An  isotropic  hose  wall  is  assumed.  This  geometry  matches 


that  of  an  experimental  test  which  will  be  described  later.  Axisyminetric 


fluid"  finite  elements  are  used  to  model  the  inner  fluid  and  axisymmetric 


thick  shell  elements  with  complex  Young's  modulus  are  used  to  model  the 
hose  wall.  The  finite  element  idealization  and  nodal  points  are  also 


shown  in  Figure  1 


The  effect  of  the  hydrophones  on  the  pressure  pattern  in  the  inner 


fluid  is  shown  in  Figure  2.  Plotted  along  the  length  of  the  array  are 
the  differences  between  pressure  (AP)  with  and  without  the  hydrophones. 
Differences  between  radial  displacements  (AR)  and  between  axial  displacements 
(AZ)  are  also  plotted.  All  curves  are  normalized  to  the  maximum  value  of 
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P,  R,  or  Z in  the  array.  The  low  level  curve  (maximum  scattering  2 
percent)  Is  for  hydrophones  with  cross  section  0.25  of  that  of  the  fluid. 
The  high  level  curve  (maximum  scattering  10  percent)  is  for  hydrophones 
with  cross  section  0.75  of  that  of  the  fluid.  As  would  be  expected,  the 
hydrophones  scatter  the  pressure  and  radial  waves,  but  not  the  longitudinal 
waves . 

Figures  3 and  4 show  the  predicted  complex  input  impedance  of  the 

7 8 

module  for  Young's  modulus  3.39x10  and  1.0x10  , respectively.  Note  that 
the  curves  are  similar  in  shape  with  the  frequency  of  the  antiresonance 
higher  for  a large  Young's  modulus.  A typical  set  of  predicted  results 
are: 

Y - 3.39xl07  (1  - j .25)  Newt/met2 
f - 100  Hz 
F^n  » 6.28  Newtons 

- (-.17,  -,63)xl0  4 met. 

UQut  - (-.58,  .25)xl0-4  met. 

■ pressure  at  hydr.  #1 

- (-.29,  -.18)xl05  ubar 

- (-.42,  -.26)  psi 

- radial  displacement  at  hydr.  #1  (bulge  wave  amplitude) 

- (.19,  .22)xl0-5  met. 
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Figure  5 Is  a frame  from  a motion  picture  showing  a time  history  of  a 
vibrating  module.  The  radial  motion  is  exaggerated  10:1  over  the  axial 
motion.  Note:  a)  the  existence  of  a bulge  wave,  b)  the  motion  of’  the 

fluid  relative  to  the  hose  wall,  and  c)  the  vibration  in  thickness  of 
the  hose  wall,  indicating  a necessity  to  use  thick  shell  elements  to  model 
the  wall.  Figure  5 shows  that  the  model  is  simulating  the  proper  phenomena 
in  axial  vibration.  In  order  to  determine  the  validity  of  the  model  in 
simulating  an  actual  towed  array  module,  the  theory  must  be  compared  with 
controlled  laboratory  experiments.  Comparison  with  sea  tests  is  of  limited 
value  because  several  noise  phenomena  occur  simultaneously  and  generally 
hydrophone  output  is  the  only  measured  quantity.  A simple  experimental 
setup  is  shown  in  Figure  6.  It  is  intended  to  use  this  experiment  to 
determine  radial  hose  wall  motion,  input  impedance  and  hydrophone  output  as 
a function  of  frequency  for  axial  and  transverse  vibrations.  The  adequacy 
of  various  finite  element  models  can  then  be  determined  by  comparison  with 
these  results. 

Some  preliminary  experimental  results  are  shown  in  Figure  7.  The 
vertical  axis  is  the  ratio  of  radial  displacement  at  a point  10  inches 
from  the  front  of  the  module  to  axial  displacement  at  the  front  of  the 
module.  The  curve  labeled  "computed"  represents  finite  element  predicted 
results.  The  curve  labeled  "experiment"  was  obtained  from  the  apparatus 
shown  in  Figure  6.  This  curve  is  given  for  a region  where  bulge  waves 
were  actually  measured  (70-200  Hz).  Outside  this  region,  the  component 
of  transverse  motion  was  comparable  to  the  bulge  wave.  Although  the  curves 
have  the  same  shapes,  the  levels  differ  by  as  much  as  50  percent.  Further 
analysis  of  both  the  computed  and  measured  results  will  be  performed  in 
an  attempt  to  discover  the  source  of  the  discrepancies. 
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FLEXURAL  DISC  MODELING 
R.W.  Dunham 

Naval  Underwater  Systems  Center 


INTRODUCTION 


Flexural  disc  transducer  elements  are  particularly  useful  in  low 
frequency  active  sonar  systems  being  capable  of  generating  relatively 
low  frequencies  while  maintaining  reasonable  physical  dimensions.  They 
also  are  used  in  line  hydrophone  arrays  and  are  amenable  to  application 
as  drivers  in  proposed  designs  of  Helmholtz  resonators.  In  the  present 
talk,  some  results  obtained  from  a special  purpose  computer  program  (see 
Reference  1)  for  modeling  axisymmetric,  multilaminate,  piezoelectric 
flexural  discs  including  the  support  structure,  or  hinge,  will  be  presented 
The  case  considered  involves  a disc  with  no  radiation  loading  (i.e.,  in 
air).  The  results  show  an  optimum  hinge  geometry  for  the  class  hinge 
studies  as  far  as  the  coupling  coefficient  is  concerned. 

THE  GENERAL  FLEXIBILITY  MODEL 


The  program  is  based  on  a general  flexibility  model  for  a multi- 
laminate disc  with  a hole  at  the  center  (solid  disc  is  a special  case). 
The  disc  is  broken  up  into  circular  contours  which  serve  as  mechanical 
"ports"  to  render  a matrix  description  of  the  disc's  elastic  behavior. 
The  program  itself  is  divided  into  three  basic  parts:  Subprogram  FLEX, 

Subprogram  BOUCO,  and  Subprogram  DYANA. 

Figure  1 shows  a disc  with  typical  circular  contours  and  the  port 
matrix  equation  relating  a vector  containing  the  generalized  port  forces 
to  a vector  containing  the  generalized  port  displacements.  In  addition 


the  vectors  contain  variables  involved  in  the  boundary  conditions  at  the 
inside  and  outside  edges  of  the  disc.  Hence,  the  matrix  [D]  is  an 


expanded  flexibility  matrix.  It  is  evaluated  by  Subprogram  FLEX 


The  general  dynamic  equations  for  the  n-port  disc  model  are 


Here  the  matrix  [D]  is  the  reduced  flexibility  matrix  after  the  boundary 
conditions  were  applied  by  Subprogram  BOUCO.  (Ft)  is  the  vector  of  total 
force  given  here  as  the  sum  of  the  vector  of  applied  generalized  force 
and  load  terms  involving  the  generalized  displacement  vector.  The  matrix 
[A]  relates  the  generalized  displacements.  Subprogram  DYANA  performs  the 
necessary  matrix  algebra  to  get  (W)  on  the  left-hand  side  only. 


Finally,  DYANA  solves  for  (W) 


ASSUMPTIONS  IN  THE  MODEL 


There  are  two  basic  assumptions  made  in  developing  the  model 
1.  The  strains  S.  and  S~  are  assumed  given  by: 


where  Z is  the  axial  coordinate,  r is  the  radial  coordinate,  0 is  the 


slope  of  the  deflection,  and  Y is  the  uniform  radial  displacement 


2.  The  shearing  strain  (if  Included)  is  assumed  to  be  parabolic  as 


MMMMV 


I 


in  simple  beam  bending. 


MODEL  FOR  A SIMPLE  DISC  WITH  HINGE 

The  model  for  a bilaminar  disc  of  PZT-4  with  a simple  flange  hinge  is 
shown  in  Figure  2.  The  disc  is  modeled  as  one  radial  section  since  it  is 
uniform  radially.  The  two  laminates  of  PZT-4  are  assumed  identical  in 
material  properties.  The  hinge,  however,  requires  two  radial  sections; 
one  for  the  thick  region  and  one  for  the  flange.  Also,  three  laminates 
are  required  for  the  hinge  as  shown.  The  disc  portion  is  broken  up  into 
twenty  circular  contours.  The  hinge  requires  only  one  contour  at  its 
inside  edge,  this  supplying  the  flexible  boundary  conditions  for  the  outside 
edge  of  the  disc  portion  when  the  disc  is  united  to  the  hinge  in  a BOUCO 
run.  The  outside  edge  of  the  flange  is  taken  as  clamped. 

Figure  3 shows  the  disc  and  hinge  with  the  geometrical  parameters 
pertaining  to  the  model.  Both  resonant  frequency  and  effective  coupling 
coefficient  were  calculated  vs.  the  relative  flange  length,  a/(a+b),  for 
various  relative  flange  thicknesses,  tM/t.  These  plots  are  shown, 
respectively,  in  Figures  4 and  5.  Note  that  for  the  very  thick  flange, 
a/ (a+b)  ■ .9450,  there  is  very  little  dependence  on  a/(a+b)  as  expected, 
the  values  remaining  close  to  those  for  a uniform,  full  thickness  hinge 
clamped  at  its  outer  edge.  It  is  interesting  to  note  the  behavior  of 
the  curves  for  very  thin  flanges.  As  the  flange  becomes  shorter  and 
shorter  the  curves  first  appear  to  be  approaching  values  for  a uniform, 
full  thickness  hinge,  simply  supported  at  its  outside  edge  (shown  as  S.S. 
on  the  graph)  but  finally  the  curves  take  off  toward  the  values  for  a 


16 
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uniform,  full  thickness  hinge  clamped  at  its  outside  edge,  (shown  as 
clamped  on  the  graph).  A little  reflection  indicates  this  is  in  agreement 
with  one's  intuition. 

CONCLUSIONS  AND  FUTURE  PLANS 

This  computer  program  proved  very  efficient  for  modeling  flexural 
discs  with  the  flexible  supports.  A very  large  number  of  cases  were 
run  to  generate  the  curves  shown  with  relatively  little  computer  time 
and  relatively  few  manhours  invested.  The  study  showed  that  an  optimum 
hinge  geometry  exists  for  the  coupling  coefficient,  i.e.,  for  tjt  - .2100 
and  a/ (a+b)  - .65.  It  is  also  interesting  that  the  resonant  frequency  is 
fairly  independent  of  the  flange  length  for  this  geometry.  It  should  be 
pointed  out  that  this  consideration  ignores  the  stress  distribution  at 
present. 

Future  plans  consist  of  the  following. 

1.  Actual  flexural  discs  with  supports  similar  to  the  above  are 
presently  being  constructed.  These  will  be  carefully  measured  generating 
experimental  curves  corresponding  to  the  ones  shown  here.  Corresponding 
computer  runs  will  be  made  to  compare  with  the  experimental  results  thus 
giving  a thorough  check  on  the  program. 

2.  One  shortcoming  of  the  present  program  is  the  inability  to 
output  the  stress  profile  throughout  the  disc  and  hinge.  This  deficiency 
will  shortly  be  rectified  when  this  capability  is  added. 

3.  The  model  is  to  be  modified  into  an  effective  two-port  model  to 
supply  an  effective  ABCD  matrix  for  use  with  Dave  Porter's  array  programs. 
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SUBPROGRAM  FLEX  CALCULATES  [o]  FOR  THE  MODEL 


GENERAL  FLEXIBILITY  MODEL 


figure  1 


USE  OF  ORTHO"  AL  FIELDS  IN  MAGNETOSTRICTIVE  SCROLL  TRANSDUCERS 

C.  L.  LeBlanc  and  C.  H.  Sherman 
Naval  Underwater  Systems  Center 

ABSTRACT 

A linear  magnetos trictive  theory  is  used  to  investigate  the 
alternating  small  signal  characteristics  of  a thin-wall  scroll  trans- 
ducer polarized  and  driven  with  orthogonal  magnetic  fields . The  concept 
of  orthogonal  excitation,  advocated  by  researchers  in  the  International 
Nickel  Co.,  is  studied  from  this  linear  viewpoint  and  values  for  the 
effective  electromechanical  coupling  coefficient  are  presented  in  terms 
of  the  material  coupling  coefficient  derived  by  applying  colinear 
polarization  and  excitation  magnetic  fields  along  the  axis  of  principal 
strain.  Although  the  effective  coupling  was  found  to  depend  quite 
heavily  on  the  orientation  of  the  resultant  magnetic  fields,  no  case  was 
encountered  where  the  effective  coupling  exceeded  k^.  However,  it  is 
expected  that  orthogonal  excitation  will  have  genuine  advantages  under 
large  signal  driving  conditions  and  two  special  cases  of  high  amplitude 
drive  will  be  discussed;  linearized  drive  and  unpolarized  drive. 

Magnetostrictive  scroll  transducers  have  been  in  use  for  over  twenty 
years.  In  all  this  time  the  first  new  approach  for  driving  these  scrolls 
is  the  orthogonal  drive  concept  invented  by  AR  Edson  at  INCO.  Investi- 
gation of  the  orthogonal  drive  concept  shows  that  the  coupling  coef- 
ficient for  an  orthogonally  driven  scroll  cannot  exceed  the  coupling  for 
a conventionally  driven  scroll  but  that  high  amplitude  orthogonal  drive 
gives  better  linearity  and  more  acoustic  power  than  high  amplitude  con- 


ventional drive. 


Figure  1 shows  a thin  wall  scroll  transducer  with  two  windings;  an 


exterior  toroidal  winding,  which  is  the  one  used  in  conventional  scrolls, 
and  an  interior  solenoidnal  winding  which  is  completely  surrounded  by 
magnetostrictive  material.  Bias  and  alternating  currents  can  be  super- 
imposed in  both  windings . Current  it  in  the  toroidal  winding  produces 
a closed  circumferential  magnetic  field  - H^.  Current  ig  in  the 
solenoidal  winding  produces  an  axial  magnetic  field  - - which  forms  a 

closed  path  in  the  material  as  outlined  by  the  connecting  arrows.  This 
latter  field  is  transverse  to  the  circumferential  field;  thus,  the 
term  "orthogonal  drive".  The  resultant  magnetic  field  is  and  is 
inclined  at  an  angle  0 to  the  field.  To  gain  a better  understanding 
of  orthogonal  drive,  the  scroll  transducer  was  analyzed  using  linear, 
small  signal  magnetostrictive  theory^.  The  objective  was  to  determine 
the  effective  electromechanical  coupling  coefficient  and  to  see  if 
orthogonal  drive  provided  higher  coupling  than  conventional  drive. 

It  is  conventional  to  write  the  linear,  small  signal  equations  of 
state  in  a coordinate  system  with  the  x^  direction  parallel  with  the  bias 
magnetization  direction  in  the  material.  In  Figure  2,  the  resultant  bias 
magnetization  vector  (H^)  is  inclined  at  an  angle  0 such  that  the  tangent 
of  0 is  the  ratio  of  the  physically  applied  orthogonal  bias  fields , H| 
to  H^.  The  first  step  in  the  analysis  is  to  write  the  constitutive 
equations  of  state  for  polarized  polycrystalline  material  in  the  unprimed 
coordinate  system.  Since  the  unprimed  system  is  rotated  by  the  angle  0 
to  the  primed  system,  which  fits  the  scroll  geometry  and  by  which  the 
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scroll  motion  can  be  described,  then  the  original  equations  of  state 
must  be  transformed  to  the  primed  system  for  practical  analysis.  The 
transformation  is  accomplished  using  tensor  analysis  and  the  transformed 
equations  which  are  needed  in  the  prime  coordinate  system  are  those 
shown  in  Figure  3. 

These  particular  equations  apply  to  a simple  magnetos trictive  ring 
transducer  with  magnetic  isotropy  and  relate  strain  S,  and  stress  T to 
magnetic  intensity  H and  magnetic  induction  B.  Assuming  that  the  mech- 
anical stress  and  strain,  as  well  as  the  magnetic  intensity  and  induction 
are  uniform  throughout  the  material,  then  the  static  and  dynamic  coupling 
coefficients  are  equal. 

The  coupling  coefficient  is  evaluated,  using  an  energy  approach 
defined  in  the  IEEE  Standards  on  Piezoelectric  Crystals,  by  calculating 
the  piezomagnetic  strain  and  induction  energy  densities.  The  specific 
energies  are  found  by  multiplying  equation  (1)  by  the  only  existing  stress 
component,  the  T^  circumferential  stress,  and  by  multiplying  equations 
(2)  and  (3)  by  their  respective  alternating  magnetic  field  components, 
and  H^.  The  results  are  shown  in  Figure  4.  Equation  (1)  gives  the 
magnitude  of  the  effective  coupling  in  terms  of  the  unprimed  coefficients, 
the  ones  usually  measured  and  published  as  material  constants.  The 
relationships  between  the  primed  and  unprimed  coefficients  are  found  by 
comparing  the  appropriate  bracketed  expressions.  In  equation  (1),  H' 
is  the  ratio  of  the  magnitudes  of  the  alternating  magnetic  field  components, 
to  H^.  If  the  following  assumptions  are  made: 
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(1)  the  material  Is  isotropic  elastically, 

(2)  there  is  no  change  in  volume  upon  magnetization,  and 

(3)  the  elastic  and  piezomagnetlc  shear  coefficients  have 
relationships  similar  to  their  piezoceramic  counterparts, 

then  equation  (1)  reduces  to  the  calculable  form  given  by  equation  (2). 

The  material  coupling  coefficient  is  associated  with  a conventionally 
driven  scroll.  Many  cases,  ranging  from  H'  equals  zero  to  H'  equals 
infinity,  are  possible  in  equation  (2);  only  specific  cases  will  be 
commented  upon. 

Figure  5 shows  the  coupling  coefficient  ratio  plotted  as  a function 
of  the  bias  direction  0 for  various  driving  signal  ratios  H'.  The  upper 
left  plot  shows  that  the  coupling  factor  varies  as  the  cosine  of  0 when 
H'  equals  zero  (or,  when  the  only  active  driving  field  is  the  field). 
For  0 equals  zero,  the  driving  field  is  parallel  with  the  strain  and 
the  bias  vector,  and  maximum  coupling  is  achieved.  This  case 
represents  the  conventionally  driven  scroll.  For  0 equals  ninety  degrees, 
the  driving  field  is  transverse  to  the  magnetic  bias  vector,  which  is  now 
in  the  xj^  direction,  and  the  result  is  zero  coupling.  The  plot  in  the 
upper  right  hand  corner  shows  that  the  coupling  factor  varies  as  one- 
half  the  sine  of  0 when  H'  equals  infinity  (or,  when  the  only  active 
driving  field  is  now  the  H|  field).  For  0 equals  zero,  the  driving  field 
is  transverse  to  the  bias  vector  and  the  result  is  again  zero  coupling. 
For  0 equals  ninety  degrees,  the  driving  field  is  parallel  with  the  bias 
vector  in  the  xj  direction,  but  is  transverse  to  the  strain, 
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normalization  criteria  (where  subscript  s designates  saturation) , is 
used  in  equation  (1)  to  derive  equation  (3).  Equation  (3)  gives  the 
useful  normalized  strain  in  terms  of  the  normalized  Cartesian  components 
of  magnetic  induction.  We  will  now  use  this  formula  in  considering  the 
case  which  provides  linear  high  amplitude  operation,  as  outlined  in 
Figure  7.  Let  the  components  of  Induction  be  given  by  equations  (1)  and 
(2),  where  bpx  and  b are  bias  components  and  u(t)  and  v(t)  are,  in 
general,  arbitrary  functions  of  time.  Inserting  these  components  in  the 
strain  formula,  given  in  Figure  6 yields  equation  (3).  Applying  the 
condition  shown  in  equation  (4)  eliminates  the  bracketed  term  with  the 
squared,  time  dependent  functions  on  the  extreme  right-hand  side  of 
equation  (3).  The  expression  for  the  strain  then  reduces  to  the  form 
given  in  equation  (5).  Since  both  bracketed  terms  in  equation  (5)  are 
constants,  the  strain  is  a linear  function  of  the  magnetic  induction 
u(t).  When  the  induction  is  a sinusoidal  function  of  time,  the  bias 
components  can  be  adjusted  to  maximize  the  alternating  peak-to-peak 
strain  to  a value  of  the  square  root  of  two  times  saturation  strain.  In 
a conventionally  biased  scroll  the  peak-to-peak  strain  cannot  exceed 
saturation  strain  and  there  is  considerable  second  harmonic  distortion. 


Figure  8 shows  the  other  case  of  Interest,  unpolarlzed  orthogonal  drive. 
Let  the  components  of  induction  be  given  by  equations  (1)  and  (2) . These 
conditions  cause  the  resultant  magnetic  induction  vector  to  be  constant 
in  magnitude  and  to  rotate  at  an  angular  frequency,  omega  (w) . Inserting 
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saturation  induction  makes  bQ  equal  unity,  and  then  the  strain  has 
a fixed  component  of  one-quarter  saturation  strain  with  an  alternating 
component  of  three-quarters  saturation  strain.  Therefore,  the  peak-to- 
peak  strain  is  three-halves  saturation  strain.  When  a conventional 
scroll  is  driven  unpolarlzed,  the  peak-to-peak  strain  cannot  exceed 
saturation  strain.  In  both  cases,  the  strain  varies  at  twice  the 
frequency  of  the  driving  induction. 

In  summary,  we  note  that  the  coupling  coefficient  for  an  orthogonally 
driven  scroll  cannot  exceed  the  coupling  for  a conventionally  driven 
scroll.  We  also  note  that  polarized  orthogonal  drive  gives  higher  strain 
amplitude  with  better  linearity  than  conventional  polarized  drive  and 
that  unpolarlzed  orthogonal  drive  gives  higher  strain  amplitude  than 
conventional  unpolarlzed  drive.  For  the  latter  cases,  orthogonal  drive 
has  the  potential  for  about  3 dB  more  acoustic  power. 
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COMPARISON  OF  TWO  METHODS  FOR  ANALYZING  FREE-FLOODED 
RING  ARRAYS  OF  TRANSDUCERS 


D.T.  Porter 

Naval  Underwater  Systems  Center 


Two  computer  programs  used  for  predicting  the  transmitting  behavior 


of  coaxial  free-flooded  arrays  of  transducer  rings  are  discussed  in  this 


paper. 


Figure  1 shows  a cross  section  of  an  array  of  three  transducer  rings, 


labelled  T^,  T^,  T^.  In  the  gaps  between  transducers,  and  at  the  ends 


of  the  array  are  placed  "water  rings",  having  the  thickness  of  the  trans- 


ducer rings  and  the  density  of  water.  Mutual  coupling  coefficients  between 


all  transducer  and  water  rings  are  calculated  in  conventional  cylindrical 


coordinates  by  Rogers'  equations  on  the  outer  and  inner  surfaces  of  an 


infinite  rigid  cylinder.  Thus,  the  inner  and  outer  surfaces  of  the 


transducer  rings  are  connected  by  the  action  of  the  water  rings. 


Unfortunately,  the  formulation  of  this  "water  ring"  method  requires  the 


extension  of  the  cylindrical  surface  past  the  end  water  rings  to  be  rigid. 


Because  of  this  mathematical  approximation,  the  predicted  pattern  will  be 


inaccurate  near  the  array's  axis. 


In  Figure  2 we  show  that  for  M transducer  rings  and  N water  rings 


array  equations  for  M+N  radiators  can  be  derived,  as  in  equation  (1) , 


F is  the  driving  force  vector,  Z is  the  mutual  impedance  matrix,  and  V 


is  the  velocity  vector.  The  velocity  vector  is  solved  for  by  inverting 


the  Z matrix,  as  in  equation  (2).  The  NUSC/NL  program  is  presently  limited 


to  M+N  less  than  or  equal  to  60.  By  breaking  the  program  into  three 
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separate  parts  this  limitation  could  be  raised  to  200.  However,  the 
largest  array  of  coaxial  free-flooded  rings  of  which  we  have  knowledge  is 
18. 

An  important  assumption  of  this  program  is  that  the  water  rings  with 
their  discrete  velocities  will  adequately  simulate  the  continuous  velocity 
variations  that  actually  exist  in  the  gaps  between  transducer  rings  and 
at  the  ends  of  the  array. 

Calculations  have  been  done  for  arrays  with  various  numbers  of  water 
rings  per  gap  and  per  end.  If  the  gap  is  much  less  than  a wavelength  and 
also  less  than  the  transducer  height,  then  one  or  two  water  rings  per  gap 
will  suffice.  However,  on  the  extension  of  the  cylindrical  surface  at  the 
ends  of  the  array  the  velocity  normal  to  the  cylindrical  surface  is  changing 
rapidly  in  magnitude  and  phase.  As  yet,  a satisfactory  method  for 
determining  the  required  number  and  size  of  end  water  rings  has  not  been 
determined . 

The  program  as  written  requires  the  inner  and  outer  velocities  of  any 
ring  to  be  equal  and  out  of  phase,  and  the  top  and  bottom  edges  of  each 
ring  have  zero  normal  velocity.  The  program  contains  lumped  circuits  for 
simple  ceramic  and  magnetostrictive  rings.  The  output  includes  source 
level,  directivity  index,  beam  patterns,  radiation  impedance  data,  currents, 
voltages,  powers,  and  efficiency. 

The  second  program  (ref.  b)  is  based  upon  the  Helmholtz  Integral 
formulation.  In  this  program  we  combined  the  NRL  Multiple  Ring  SHIP 
Program  (courtesy  of  Peter  Rogers)  with  the  lumped  element  transducer 
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model  of  the  Water-Ring  Program.  The  computer  input  and  output  of  the 
two  programs  (i.e.,  Rogers'  and  Porter's)  contain  nearly  the  same  items. 

In  the  SHIP  model,  the  transducer  rings  are  subdivided  into  bands 
on  their  four  surfaces,  as  is  shown  in  Figure  3.  The  SHIP  model  calculates 
mutual  impedance  coefficients  between  the  iC^  and  rings  of  an  array  by 
imposing  a unit  velocity  on  the  it^1  ring,  and  calculating  the  pressure  on 
the  jth  ring.  For  any  ring,  the  ratios  of  the  inner,  outer,  and  edge 
velocities  must  remain  fixed  over  the  array.  The  outer  velocity  is 
treated  as  the  reference  velocity.  These  reference  velocities  are  determined 
by  the  matrix  - inversion  solution  of  the  force- impedance  - velocity 
equations  shown  in  Figure  2.  In  calculating  mutual  impedance  coefficients 
(Z^),  the  Helmholtz  Integral  model  takes  into  account  the  scattering 
produced  by  all  the  rings.  If  each  of  the  N transducers  have  M bands,  then 
NxM  equations  are  required  to  solve  for  the  pressures  on  the  NxM  bands. 

Note  that  in  the  Water-Ring  Program  the  calculation  of  Z^  between  bands 
on  an  infinite  cylinder  involved  only  bands  i and  j. 

We  now  examine  the  matrix  of  a four  ring  array,  using  the 
Helmholtz  Integral  Program.  The  rings  have  a mean  radius  of  2.58  inches, 
a height  of  0.75  inch,  a thickness  of  0.25  inch,  and  a small  gap  between 
rings  of  0.06  inch.  The  array  is  resonant  in  water  around  10  kHz;  the 
Z^  presented  are  for  4 kHz.  The  inner  and  outer  velocities  on  the  ring 
surface  are  assumed  to  be  equal  and  out  of  phase;  the  top  and  bottom  edge 
velocities  were  assumed  to  be  zero. 

Figure  4 shows  the  1-1  and  1-2  terms  of  the  Z matrix,  for  various 
numbers  of  side  and  top  bands.  Note  that  as  the  number  of  bands  is 
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increased  the  convergence  of  the  matrix  elements  is  slow,  particularly 
for  the  resistive  elements,  and  R^*  A free-flooding  ring  transducer 
is  basically  a dipole  radiator,  with  a large  amount  of  cancellation  of 
the  radiation  from  the  two  sides.  Many  bands  may  be  needed  to  predict  this 
cancellation.  Although  the  ring  thickness  at  4 kHz  is  only  about  A/60, 
several  bands  are  needed  on  the  top  and  bottom  edges,  because  the  pressure 
is  changing  very  rapidly  across  these  surfaces.  Note  the  and  for 

I (top  bands)=2  and  J(side  bands)=6.  For  1=2  and  J=8  there  is  only  about 
a two  percent  decrease  in  R^  and  R^  compared  to  1=2  and  J=6.  However, 
for  the  1=4,  J=6  case  there  was  about  a nine  percent  decrease  in  R^  and 
R^2  compared  to  1=2,  J=6.  Also  we  note  that,  for  all  the  cases  in  Figure 
4,  R^2  > R^,  so  that  the  mutual  resistance  between  rings  1 and  2 was 
actually  greater  than  the  self-resistance  of  ring  1.  In  conventional 
baffled  arrays,  it  is  unheard  of  for  the  mutual  resistance  between  two 
radiators  to  exceed  the  self-resistance  of  one  of  them.  Apparently,  the 
pressure  cancellation  in  the  R^ , calculation  was  less  than  for  U | f . At 
this  frequency,  > X^.  and  >:*  R^  or  Rj^.  However,  if  the  ring's 

Thevenin  internal  reactance  nearly  cancelled  the  radiation  reactance,  then 
the  condition  of  R^  > R-q  c°uld  cause  the  array's  velocity  control  to  be 
erratic. 

The  multiple  ring  SHIP  program  operating  on  the  NUSC  Univac  1108 
can  use  up  to  96  total  bands  in  an  array.  For  1=4,  J=8,  there  are  24  bands 
per  ring,  which  would  limit  the  maximum  number  of  rings  to  only  4!  Using 
smaller  I and  J enables  one  to  handle  larger  arrays,  but  with  less  accuracy. 
However,  there  seems  no  reason  why  the  required  I and  J should  get  larger 
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as  the  gap  width  becomes  larger.  Recall  that  with  the  Water-Ring  Program 
the  calculation  becomes  more  difficult  as  gap  width  increases. 

Figure  5 compares  calculated  and  measured  (ref.  c)  transmitting  current 
responses  for  this  same  four  ring  array.  At  most  frequencies,  the 
Helmholtz  Integral  program  agrees  more  closely  with  the  measurement  than 
does  the  Water-Ring  program.  The  gaps  between  rings  are  about  A/100  at 
10  kHz,  and  the  gaps  are  also  much  less  than  the  ring  height,  so  that  the 
water  rings  should  well  simulate  the  action  in  the  gaps.  However  the 
whole  array  is  only  slightly  longer  than  A/2  at  10  kHz.  Therefore,  the 
array  end  effects  are  important,  and  the  Water-Ring  Program  is  not 
reliable  in  modeling  end  effects. 

The  calculations  did  not  take  into  account  any  of  the  array's 
mounting  structure.  Slightly  off  of  each  end  of  the  array  was  a fiberglass 
ring  which  supported  the  wires  that  excited  the  magnetostrictive  material. 
The  fiberglass  rings. must  have  had  some  effects  on  radiation  around  tlie 
ends  of  the  array,  and  are  partly  responsible  for  the  difference  between 
the  measurement  and  the  Helmholtz  Integral  calculations.  Other  causes 
for  this  difference  are  inaccuracies  in  the  transducer  circuit  model,  in 
the  Z matrix,  and  in  the  assignment  of  the  relative  velocities  on  the 
different  ring  surfaces. 

Figure  6 compares  calculated  for  the  two  mathematical  models  for 

the  same  four  ring  array.  The  calculated  for  the  Water-Ring  model  is 
lower  for  all  frequencies,  especially  10-14  kHz.  This  difference  is 
calculated  is  consistent  with  the  calculated  beam  patterns. 
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Figures  7,  8,  9,  and  10  compare  the  measured  patterns  with  those 
calculated  at  6.0,  9.5,  12.0,  and  14.0  kHz.  For  most  of  these  frequencies 
all  three  of  these  patterns  agree  well  except  near  the  array  axis.  In 
general,  the  Helmholtz  Integral  patterns  are  more  directive  than  the 
measured  patterns,  which  are  in  turn  more  directive  than  the  Water  Ring 
patterns. 

One  of  the  objectives  of  this  project  is  to  establish  which  types  of 
arrays  are  best  suited  to  each  model.  Two  parameters  are  used  to  categorize 
an  array.  One  is  simply  the  number  of  rings.  The  second  is  the  ratio  of 
the  gap  width  to  the  ring  height. 

Figure  11  shows  some  of  our  conclusions  with  respect  to  the  usefulness 
of  the  two  models.  This  chart  is  based  on  some  intuition  and  a limited 
amount  of  array  calculations  (not  presented  here) . The  chart  suggests  that 
the  Water-Ring  program  gives  poor  results  for  sparse  arrays,  and  the 
Helmholtz  Integral  program  deteriorates  for  arrays  of  more  than  10  rings. 

In  those  cases  where  either  program  is  useable,  the  Helmholtz  Integral 
program  should  be  more  accurate.  Presently  neither  program  can  accurately 
predict  the  behavior  for  a large,  sparse  array.  Hopefully  this  situation 
can  be  improved  by  new  approximations,  probably  based  on  the  Helmholtz 
Integral  formulation. 
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WATER  RING  MODEL  FOR  3-RING  ARRAY 


EACH  RING  IS  ACOUSTICALLY  SUBDIVIDED  WITH 
3 TOP/BOTTOM  BANDS  AND  4 INSIDE/OUTSIDE 
BANDS  FOR  ANALYSIS  USING  HELMHOLTZ 
INTEGRAL  PROGRAM. 


CROSS  SECTION  OF  2-RING  ARRAY 


FIGURE  3 
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TRANSMIT  CURRENT  RESPONSE  FOR  A 4- RING  ARRAY 


DIRECTIVITY  INDEX  FOR  A 4-RIN6  ARRAY 
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AREAS  OF  USEFULNESS  OF 
HELMHOLTZ  INTEGRAL  AND  WATER 
RING  PROGRAMS 
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MECHANICAL  ACOUSTIC  RADIATOR 
M.L.  Rumerman 

Naval  Ship  Research  and  Development  Center 

This  paper  outlines  the  analysis  of  a mechanically-driven  acoustic 

radiator.  The  ultimate  goal  is  to  provide  a guide  by  which  the  device 

can  be  designed  to  operate  broadband  over  a specified  frequency  range. 

The  radiator  consists  of  a thin  cylindrical  shell  and  the  mechanical 

excitation  is  generated  by  diametrically-opposed  flatted  rollers 

(Figure  1),  which  are  made  to  traverse  the  circumference  of  the  inner 

surface  of  the  shell  at  constant  angular  speed.  When  a flat  comes  down 

upon  the  shell,  an  impulsive  force  is  delivered,  and  the  vibration  thus 

generated  results  in  sound  radiated  into  the  acoustic  medium. 

Were  the  rollers  to  move  without  slipping  on  the  shell,  only 

discrete  frequencies  would  be  generated,  due  to  the  kinematic  regularity. 

Slipping  results  in  a continuous  spectrum.  Let  A be  the  average  interval 

between  successive  impacts.  Then  A = A + £ where  A is  the  deterministic 

time  between  impacts  (dependent  upon  geometry  and  speed)  and  c,  is  the 

average  slip  time  between  successive  impacts.  It  can  be  shown  that  when 

t • where  fmin  is  the  lowest  frequency  of  interest,  then  the  impacts 

min 

can  be  considered  uncorrelated  and  the  resulting  spectral  density  of 
radiated  power,  S(co),  can  be  given  by 

S (oo)  = co  > 2nf 

T min 
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where  F(co)  is  the  average  spectral  density  of  total  energy  radiated  by  a 
single  impact. 

F(w)  can  be  calculated  by  a modal  analysis  of  the  shell  vibration 

and  radiation.  It  is  found  that  the  spectral  density  of  the  excitation 

will  be  broadband  when  2tt  f e < 1,  where  f is  the  highest  frequency 

max  max  J 

of  interest  and  e is  the  duration  of  an  impact.  The  radiated  power 
spectrum  will  then  be  broadband  if  the  shell  can  be  made  to  respond 
uniformly  over  the  frequency  range. 

Figure  2 represents  idealized  modal  transfer  functions  - from 
excitation  spectrum  to  radiated  power  spectrum  - as  a function  of  frequency. 
Assuming  that  the  peaks  are  equal  and  that  adjacent  curves  cross  at  their 
mutual  3 dB  down  points  and  that  adjacent  modes  are  independent,  it  can  be 
shown  that  the  effect  is  to  give  a broadband  response  over  the  frequency 
range,  as  indicated  by  the  dashed  line.  While  it  is  impossible  to  design 
a shell  having  the  acoustical  features  described  above,  something  approaching 
it  should  be  feasible  and  would  provide  efficient  broadband  response. 
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TWO  METHODS  FOR  PREDICTING 
ARRAY  DIRECTIVITY  INDEX 


D.  T.  Porter 

Naval  Underwater  Systems  Center 


ABSTRACT 


The  first  method  is  referred  to  as  the  pattern  method.  Pertinent 
equations  are: 


N 


DI 


10  log  in  M.R.A.\ 

10  y Average  Intensity  f 


(1) 


and 


fo"K 


Average  Intensity  = I 


AV 


oI(6.<i))r  sin  0d Ud <t> 


4Trr 


(2) 


Unfortunately,  to  accurately  evaluate  the  integral  of  equation  (2),  a 
great  many  pattern  points  may  be  needed. 

The  second  method  is  referred  to  as  the  radiated  power  method.  If 
the  array  under  consideration  is  a receiving  array,  we  must  mathematically 
concoct  an  analogous  and  reciprocal  transmitting  array.  Pertinent 
equations  are 


Source  Level  = 71.6  + N^  + 10  log  (Acoustic  Power) 


O) 


and 


Acoustic  Power  ' R.V 


j-1 


J j 


(4) 


We  may  find  N^^  from  (3)  if  we  know  the  source  level  on  the  MRA,  and  the 
acoustic  power  radiated  from  each  transducer.  Knowledge  of  the  source 
level  requires  the  array's  acoustic  response  at  only  one  point,  instead 
of  many  as  in  the  pattern  method.  However,  knowledge  of  the  acoustic 
power  requires  solution  of  the  array's  coupled  force  - velocity  - mutual 
coupling  equations  in  order  to  get  the  radiation  resistances  and  velocities 
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needed  in  equation  (4) . This  in  turn  requires  expressions  for  mutual 
coupling  coefficients  (Z^).  The  programming  and  debugging  of  these 
expressions  may  be  very  troublesome. 

Provided  the  Z can  be  obtained,  the  radiated  power  method  gives 
exact  results  for  N^.  An  example  was  shown  where  the  computer  time 
(Univac  1108)  for  a 1000  - element  array's  was  100  minutes  for  the 
pattern  method  and  22  minutes  for  the  radiated  power  method. 
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COMPUTER  TIME  FOR  UNIVAC  1108 
TO  CALCULATE  VIA  PATTERNS 


I 


I 


kJ 


N N 
-E 


(min)  36,000 


For  2-D  Pattern,  N = N,  *N 
P V 4 


NflN  N 

T = — ....9  c 

(min)  36,000 


For  NQ  = 180,  = 20,  and  Nt  = 1000 


T(min)  = 100  rainutes 


(Same  array  required  22  minutes  using 
Radiated  Power  Method) 
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PATTERN  METHOD  OF  OBTAINING  N 


RADIATED  POWER  METHOD 
OF  OBTAINING  N_ 


L = 71.6  + N_t  + 10  log  P 
s DI  8 ACOUSTIC 


NDI  = Ls  "71*6  " 10  log  P 
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A FINITE  ELEMENT  APPROACH  TO  ACOUSTIC 
RADIATION  FROM  ELASTIC  STRUCTURES 

M.  Knittel 
J.  Hunt 
D.  Barach 

Naval  Undersea  Center 

Past  mathematical  models  of  vibrating  structures  immersed  in  an 
infinite  fluid  medium  have  not  been  entirely  successful  because  of  the 
difficulties  involved  in  combining  the  vibrational  aspects  of  the  structure 
with  the  acoustic  radiation  properties  of  the  medium.  In  cases  where 
this  interaction  has  been  successfully  treated,  such  as  in  a combination 
of  a finite  element  program  (MARTSAM)  with  a Helmholtz  Integral  radiation 
program  (CHIEF),  the  necessity  for  evaluating  numerous  and  expensive 
integrals  limits  the  usefulness  of  these  methods.  Near  and  far  field 
pressures  are  calculated  at  a minimum  number  of  field  points  for  a limited 
number  of  frequencies.  The  approach  that  we  have  taken  to  acoustic 
radiation  from  elastic  structures  avoids  this  difficulty;  calculation  of 
the  near  and  far  fields  for  numerous  frequencies  is  an  inexpensive  process 
and  thus  allows  the  calculation  of  the  pressure  at  enough  points  to  draw 
contour  plots  of  the  pressure  field. 

The  problem  chosen  for  analysis  with  this  method  was  the  behavior  of 
a tangentially  poled,  piezoelectric,  ceramic  cylinder  immersed  in  an 
infinite  fluid  medium  (Figure  1).  This  free-flooded  cylinder  was  chosen 
because  it  is  a simple,  axisymmetric  body  in  which  circular  harmonics 
other  than  m=0  can  easily  be  electrically  excited;  that  is  a voltage  that 
varies  as  cosine  m0  can  be  applied  around  the  cylinder.  These  non-zero 
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harmonics  are  important  in  scattering  from  a single  cylinder  because  all 
circular  harmonics  will  be  excited  when  the  incident  wave  does  not 
propagate  down  the  symmetry  axis.  The  non-zero  harmonics  are  also 
important  in  the  design  of  non-coaxial  arrays  of  cylinders  since  interactions 
between  them  will  excite  non-zero  harmonics  of  the  individual  cylinders. 

Finite  element  techniques  were  used  to  model  one-half  of  a double 
stave  of  the  cylinder  and  the  appropriate  boundary  conditions  applied 
(Figure  2).  Each  nodal  point  has  three  displacement  degrees  of  freedom 
and  an  electrical  degree  of  freedom.  The  coordinate  axes  of  the  structure 
were  transformed  from  rectangular  to  a cylindrical  coordinate  system  and 
the  degrees  of  freedom  condensed  down  to  an  axisymmetric  representation. 
Cosine  m 0 voltage  dependence  was  applied  around  the  cylinder,  and  cos  m0 
and  sine  m 0 displacement  dependence  was  assumed.  Below  20  kHz,  our 
calculations  have  shown  that  this  reduction  is  valid. 

The  fluid  immediately  surrounding  the  cylinder  was  modeled  using 
finite  element  techniques  out  to  a sphere  that  just  encloses  the  structure. 
This  was  done  because  for  spherical  surfaces,  the  acoustic  radiation  Green's 
function  for  steady-state  problems  separates  when  expanded  in  spherical 
harmonics.  A closed  form  solution  for  the  radiation  impedance  that  is 
associated  with  each  spherical  harmonic  can  be  found  that  consists  of 
merely  a ratio  of  Hankel  functions.  The  finite  element  fluid  sphere 
surface  forces,  pressures,  and  normal  velocities  can  be  expanded  in 
approximations  to  spherical  harmonics  through  a simple  transformation  to 
a finite  set  of  spherical  harmonic  base  vectors.  The  spherical  harmonic 
components  of  the  pressure  at  any  point  outside  the  finite  element  fluid 
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sphere  can  be  calculated  in  terms  of  the  sphere  surface  pressure  components 
and  ratios  of  two  Hankel  functions. 

Forty-nine  acoustic  fluid  elements  were  used  to  model  one-quarter  of 
the  fluid  sphere  in  the  x-y  plane  (Figure  3) ; these  axisymmetric  elements 
are  triangular  in  cross  section  with  six  nodes  per  element.  Each  node  has 
one  degree  of  freedom,  the  pressure.  A connection  or  comparability  matrix 
was  formed  that  relates  the  normal  displacements  of  the  structure  to  the 
pressures  in  the  fluid  on  the  common  surfaces.  This  gives  a total  of 
about  150  degrees  of  freedom.  However,  by  solving  the  interior  problems 
of  the  free  vibrations  of  the  combined  structure  - fluid  sphere  with  zero 
sphere  surface  pressure,  the  number  of  degrees  of  freedom  can  be  reduced 
down  to  about  twenty-five.  In  addition,  solving  the  equations  in  this 
manner  avoids  singularities  that  occur  for  frequencies  near  eigenf requencies 
of  the  interior  problem.  Therefore,  once  the  finite  element  mass  and 
stiffness  matrices  are  formed  and  assembled  and  the  interior  problem  solved 
once,  only  matrices  of  the  order  of  25x25  need  to  be  solved  for  each 
driving  frequency.  No  acoustic  radiation  integrals  are  numerically 
evaluated,  and  the  pressures  at  a large  number  of  field  points  for  numerous 
frequencies  can  be  inexpensively  calculated. 

Figure  4.  m=0  complex  electrical  impedance.  The  real  part  of  the 

impedance  is  proportional  to  the  acoustic  power  radiated  into  the  water. 
Choosing  a frequency  of  3,000  Hz. 

Figure  5.  Structural  displacements  show  breathing  mode. 

Figure  6.  Vertical  directivity  pattern. 

Figure  7.  Near  field  pressure  contour  plot  out  to  0.4  meters  at  one 
time  step. 


67 


Figure  8.  Contour  plot  extending  out  into  the  far  field  to  3.6  meters. 


Far  field  begins  at  about  2.5  meters. 

Figure  9.  m=l  complex  electrical  impedance.  In  performing  the 
experiments,  it  was  found  that  f iberglassing  the  cylinder  lowered  the 
resonant  frequencies.  Difficulties  were  also  encountered  in  measuring  the 
magnitude  and  phase  of  the  electrical  impedance  accurately  while  driving 
the  cylinder  in  the  non-zero  circular  harmonics.  In  addition,  the  material 
parameters  used  in  our  model  are  loss-less  parameters  whose  values  in  some 
cases  are  probably  not  too  accurate.  All  of  this  contributes  to  the 
discrepancies  between  theory  and  experiment.  Look  at  4,500  Hz. 

Figure  10.  Structural  displacements  show  cosine  10  dependence. 

Figure  11.  Vertical  directivity  pattern. 

Figure  12.  Horizontal  directivity  pattern  showing  two  nodes. 

Figure  13.  Near  field  pressure  contour  plot  in  0=0  plane  showing 
central  nodal  line. 

Figure  14.  Pressure  contour  plot  extending  out  into  the  far  field. 
Figure  15.  m=2  complex  electrical  impedance. 

Figure  16.  Structural  displacements  at  8,750  Hz. 

Figure  17.  Vertical  directivity  pattern. 

Figure  18.  Horizontal  directivity  pattern  showing  four  nodes. 

Figure  19.  Near  field  pressure  contour  plot  in  0=0  plane  again 
showing  central  nodal  line. 

Figure  20.  Pressure  contour  plot  extended  out  into  the  far  field. 
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The  next  application  for  which  we  hope  to  use  this  method  is  the 
receive  or  scattering  problem  using  a free-flooded  cylinder  or  some  other 
structure.  We  also  hope  that  this  formulation  will  offer  a feasible 
method  to  solve  the  combined  equations  of  motion  in  the  time  domain.  In 
addition,  we  will  probably  use  this  method  to  design  a free-flooded 
cylinder  for  optimum  performance.  We  are  also  looking  at  new  and  more 
informative  ways  to  display  the  generated  data,  perhaps  such  as  a movie 
of  lines  or  surfaces  of  constant  pressure  that  change  with  time  and  use 
continuously  varying  colors  to  represent  magnitudes  of  pressure.  In 
addition,  we  are  trying  display  methods  that  will  show  both  the  structural 
displacements  and  the  pressure  field  and  how  they  vary  with  time  in  one 
changing  display. 
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RADIATED  PRESSURE  PATTERN  FROM  A 
TANGENTIALLY  POLED  CYLINDER 
(NEAR  FIELD) 


CIHCULAR  HARMONIC  = 0 
FREQUENCY  = 3000  Hz 
WT  = v/2 

PRESSURE  UNIT  = 1 NEWTON/(METER)2 


FIGURE  7 
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RADIATED  PRESSURE  PATTERN  FROM  A 
TANGENTIALLY  POLED  CYLINDER 
(FAR  FIELD) 


CIRCULAR  HARMONIC  = 0 
FREQUENCY  - 3000  Hz 
WT  = rr/2 

PRESSURE  UNIT  = 1 NEWTON/(METER)2 


FIGURE  8 
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FIGURE  10 


HORIZONTAL  DIRECTIVITY  PATTERN 


CIRCULAR  HARMONIC  = 1 
FREQUENCY  = 4500  Hz 

WT  = tt/2 

PRESSURE  UNIT  = 1 NEWTON/(METER) 


FIGURE  13 


RADIATED  PRESSURE  PATTERN  FROM  A 
TANGENTIALLY  POLED  CYLINDER 
(NEAR  FIELD) 


RADIATED  PRESSURE  PATTERN  FROM  A 
TANGENTIALLY  POLED  CYLINDER 
(FAR  FIELD) 


CIRCULAR  HARMONIC  = 1 
FREQUENCY  = 4S00  Hz 

WT  = rr/2 

PRESSURE  UNIT  - 1 NEWTON/(METER) 


FIGURE  14 
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RADIATED  PRESSURE  PATTERN  FROM  A 
TANGENTIALLY  POLED  CYLINDER 
(NEAR  FIELD) 


CIRCULAR  HARMONIC  = 2 
FREQUENCY  = 8750  Hz 
VVT  = ir/2 

PRESSURE  UNIT  = 1 NEWTON/(METER) 


FIGURE  19 


RADIATED  PRESSURE  PATTERN  FROM  A 
TANGENTIALLY  POLED  CYLINDER 
(FAR  FIELD) 


CIRCULAR  HARMONIC  = 2 
FREQUENCY  = 8750  H* 

WT  - r/2 

PRESSURE  UNIT  = 1 NEWTONA'MCTER) 


FIGURE  20 


APPLICATIONS  OF  THE  FINITE  ELEMENT  METHOD 
AND  COMPUTER  CRAPHICS  TO  A TRANSDUCER  ELEMENT 


J.T.  Hunt,  R,R.  Smith,  D.  Barach,  L.  McCleary 
Naval  Undersea  Center 

„ SUMMARY 

Although  many  mathematical  models  have  been  developed  for  the 
longitudinal  vibration  type  element  discussed  in  this  report,  they  have 
only  partially  succeeded  in  predicting  its  free  vibrational  response. 

The  difficulty  has  been  the  large  flexural  interaction  between  the  non- 
axisymmetric  radiating  head  and  the  electrically  driven  axisymmetric  tube, 
which  results  in  a strong  coupling  between  the  circular  harmonics  of  the 
piezoelectric  tube.  The  inability  to  predict  this  effect  has  been  the 
primary  deficiency  of  earlier  mathematical  models.  General  Dynamics/ 

Electric  Boat  Division,  in  close  cooperation  with  NUC,  has  used  the  finite- 
element  method  to  develop  a three-dimensional  model  in  which  this  deficiency 
has  been  corrected.  This  mathematical  model  is  representative  of  the  most 
advanced  techniques  which  describe  the  vibratory  characteristics  of  an 
electromechanical  sonar  transducer.  In  this  paper,  the  mathematical 
model  is  used  to  predict  the  deformation  of  the  composite  transducer  in 
;i»"- terms  of  the  free  vibrational  modes  of  its  component  parts.  Two  types  of 
computer-generated  displays  - contour  plots  of  the  predicted  deformations 
and  three-dimensional  continuously  shaded  plots  of  these  predicted  deforma- 
tions - are  used  to  graphically  illustrate  the  physical  meaning  of  the  model's 
numerical  predictions.  The  contour  plots  are  arranged  so  that  they  can 
be  directly  compared  with  the  deformations  experimentally  measured  by 
holographic  interferometry.  Predicted  deformations  are  found  to  agree 
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closely  with  those  measured  experimentally. 

The  geometrical  symmetry  of  the  transducer  consists  of  two  reflection 
planes  and  a 180-deg  rotation  axis.  By  exploiting  this  geometrical 
symmetry,  one  can  show  that  by  modeling  only  one  quarter  of  the  transducer, 
i.e.,  that  portion  of  the  transducer  situated  between  the  two  symmetry 
planes,  and  by  applying  the  correct  boundary  conditions  to  nodal  points 
located  in  the  symmetry  planes  the  correct  vibrational  response  of  the 
system  can  be  obtained. 

Although  exploitation  of:  the  symmetry  decreases  the  size  of  the 
numerical  problem  by  axmost  a factor  of  four,  the  resulting  eigenvalue 
problem  still  has  1135  displacement  degrees  of  freedom.  Hence,  each 
eigenvector  (or  the  deformation  of  the  one-quarter  model  of  the  transducer 
at  resonance)  results  in  1135  pieces  of  data,  while  the  deformation  of  the 
entire  transducer  is  represented  by  nearly  4000  numbers.  Therefore,  an 
understanding  of  the  physical  meaning  of  this  large  block  of  computer- 
generated data  can  be  achieved  only  through  its  graphic  display.  We 
have  used  two  different  types  of  computer-generated  data  displays  to 
depict  the  eigenvectors  of  the  system. 
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SELF  AND  MUTUAL  ACOUSTIC  RADIATION  IMPEDANCES 
FOR  TWO  COPLANAR  UNBAFFLED  DISKS 

A.L.  Van  Buren  and  B.J.  King 
Naval  Research  Laboratory 

SUMMARY 

The  self  ajnd  mutual  acoustic  radiation  impedances  for  two  coplanar 
unbaffled  disks  are  calculated  using  an  eigenfunction  expansion  in  terms 
of  oblate  spheroidal  wave  functions.  Terms  representing  outgoing  waves 
from  both  disks  are  included.  One  disk  is  assumed  to  be  vibrating  with 
a rotationally  symmetrical  normal  velocity  distribution,  and  the  other 
disk  is  assumed  to  be  stationary.  The  determination  of  the  expansion 
coefficients  from  these  boundary  conditions  is  facilitated  by  the  use 
of  an  addition  theorem  which  expresses  spheroidal  wave  functions  with 
reference  to  one  coordinate  frame  in  terms  of  spheroidal  wave  functions 
with  reference  to  a second  coordinate  frame.  Numerical  results  for  the 
special  cases  where  the  vibrating  disk  is  either  oscillating  uniformly 
or  vibrating  uniformly  on  only  one  side  are  presented  and  discussed. 
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INTRODUCTION 

The  surface  of  a disk  is  the  natural  coordinate  surface  E,  - 0 in 
oblate  spheroidal  coordinates  (£,n,$).  Thus  the  acoustic  radiation  from 
a disk  vibrating  harmonically  with  a known  normal  velocity-  distribution 
can  be  calculated  by  use  of  an  eigenfunction  expansion  in  terms  of  oblate 
spheroidal  wave  functions.  The  calculation  can  be  extended  to  include  a 
second,  coplanar,  vibrating  disk  by  adding  an  expansion  of  spheroidal 
wave  functions  about  the  second  disk.  An  addition  theorem  [l]  which 
expresses  spheroidal  wave  functions  with  reference  to  one  coordinate 
frame  in  terms  of  spheroidal  wave  functions  with  reference  to  a second 
coordinate  frame  can  be  used  to  facilitate  application  of  the  boundary 
conditions  on  each  disk. 

By  linear  superposition,  the  acoustic  radiation  field  of  the  two 
vibrating  disks  is  equal  to  the  sum  of  the  field  produced  when  the  first 
disk  is  vibrating  and  the  second  disk  is  stationary  and  the  field  produced 
when  the  second  disk  is  vibrating  and  the  first  disk  is  stationary.  Thus 
no  generality  is  lost  by  restricting  the  analysis  to  the  case  where  the 
second  disk  is  stationary. 

This  paper  is  concerned  with  both  the  self  acoustic  radiation  impedance 
of  the  vibrating  disk  and  the  mutual  acoustic  radiation  impedance  between 
the  vibrating  disk  and  the  stationary  disk.  The  formulas  necessary  to 
evaluate  the  self  and  the  mutual  acoustic  radiation  impedances  are  developed 
in  Sec.  2.  It  is  assumed  that  the  normal  velocity  distribution  of  the 
vibrating  disk  is  rotationally  symmetrical.  In  Sec.  3,  numerical  results 
calculated  for  the  special  cases  where  the  first  diBk  is  either  oscillat- 
ing uniformly  or  vibrating  uniformly  on  only  one  side  are  presented  and 
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discussed.  Previous  results  of  Crane  [2]  and  Hanish  [3]  vhich  were 
obtained  by  approximate  methods  that  do  not  involve  spheroidal  wave 
functions  are  also  discussed. 


ANALYSIS 


The  oblate  spheroidal  coordinates  U,n,+)  are  related  to  Cartesian 


where  d is  the  interfocal  distance,  and  where  0.*£« 


lsnsl  , and 


0*4<2*  . Here  the  surface  of  constant  £ is  an  oblate  spheroid  with  a 

2 1/2 

major  axis  of  length  U + l)  d and  a minor  axis  of  length  Cd.  The 


spheroid  for  £ ■ 0 is  a disk  of  radius  a * d/2.  The  surface  of  constant 


n is  a hyperboloid  of  one  sheet  which  intersects  the  £ » 0 surface  ii  a 

almalm  Af  4 a/l  ,1/2 


Consider  two  identical  disks  D.  and  D_  which  are  coplan ar  and  share 


a common  x axis,  as  shown  in  Figure  1.  The  radius  of  each  disk  is  a,  and 


the  center  to  center  separation  is  r^.  Spheroidal  coordinate  systems 

Cl(*l,nl’*l)  and  C2^2,n2*^2^  e8'taW-i*hed  which  contain  D1  (c^  • 0) 
and  D.  U -0),  respectively,  as  natural  coordinate  surfaces.  The 


region  exterior  to  the  disks  is  assumed  to  be  filled  with  an  infinite 
hoangeneous  fluid  of  mass  density  p and  sound  speed  c.  Disk  is 
assumed  to  be  vibrating  with  a rotationally  symmetrical  normal  velocity 
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distribution  whose  spatial  dependance  is  given  by  v(n^),  and  disk  Dg 
is  assumed  to  be  stationary.  The  time  dependence  is  assumed  to  be  har- 
monic with  the  factor  e*‘wt  being  suppressed. 

The  acoustic  pressure  produced  at  a point  r on  or  exterior  to  both 
and  Dg  can  be  expressed  as  a sum  of  outgoing  spheroidal  vaves  originat- 
ing from  both  and  Dgi 


p(r>  * X Jo 

.<M, 


(2) 


where  the  spheroidal  wave  functions  are  given  by 

- R^)(h,C)S^^(h,n)cosm^  - iR^fh.Os^h.nJcosm* 

■ ♦^)(h;e,n,4)  - l+^toC.n.*)  (3) 

- ^(h.OS^h.rOcosm*  . 

Here  the  acoustic  size  parameter  h ■ kd/2  ■ ka,  where  k * w/c  is  the  wave- 
number,  and  and  are  the  spheroidal  coordinates  of  r 

in  and  Cg.  The  spheroidal  radial  wave  functions  of  the  first  and  second 
kinds  R^^  and  R^p  and  the  apheroidal  angle  wave  functions  of  the  first 
kind  used  in  this  report  are  defined  by  Flaoswr  (k).  The  radial 
functions  R^^  and  R^  correspond  to  spherical  Bessel  functions  of  the 
first  and  second  kinds  and  are  normalized  so  that  f represents  out- 
going waves  for  time  dependence.  The  angle  functions  ^ can  be 

ft* 
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expanded  in  a series  in  the  corresponding  associated  Legendre  functions: 
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s^)Ox,n)  - f d (h|«i)p^  (n)  . (k) 

n»0,l 

where  the  prime  indicates  that  the  summation  is  over  even  or  odd  values  of 
n according  to  whether  l - m is  even  or  odd.  Terms  involving  sinm$  have 
been  omitted  from  (2)  since  the  acoustic  pressure  must  be  symmetrical 
about  both  ^ ■ 0 and  ■ 0 . 

The  coefficients  A t and  are  determined  by  the  Neumann  boundary 
conditions  imposed  on  and  : 

-ap/anj^  ■ iwpvtf^)  , on  Dx  (5) 

ap/Jng  ■ 0 , on  Dg  (6) 

where  in  spheroidal  coordinates 

(3p/3n).£  _ 0 - (l/lnla)3p/3£  . (7) 

Application  of  the  boundary  conditions  on  or  Dg  is  difficult  if  p(r) 
is  not  stressed  entirely  in  or  Cg,  respectively.  This  is  accomplished 
by  use  of  a spheroidal  addition  theorem  which  expresses  a spheroidal  wave 
function  with  reference  to  one  coordinate  frame  in  terms  of  spheroidal 
wave  functions  with  reference  to  a second  coordinate  frame.  The  general 
addition  theorem  for  spheroidal  wave  functions  which  is  derived  in  [l] 
nllowe  for  arbitrary  relative  position  and  orientation  of  the  two 
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E 
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coordinate  frames.  For  the  present  application,  the  addition  theorem 
can  be  expressed  in  the  following  form: 


where 


5 W r-  l\ i*r'‘Whl-*> 

ir,  r*»tm+l 


X l [i*”n/(2s+l)]d>-u(h|un)  £ D(s,t,r,u,m)h£2*(kr12) , 


S"U,U+1 


t»  I r-s  I 


with  the  prime  sign  indicating  that  the  summation  is  in  steps  of  two, 

where  u starts  from  0 or  1 depending  on  whether  l-m-n  is  even  or  odd, 

where  r starts  from  m or  m ♦ 1 depending  on  whether  i - m is  even  or  odd, 

and  where  s starts  from  u or  u +1  depending  on  whether  n - u is  even  or 

(2) 

odd.  The  function  h^  ' (kr1 p ) is  the  spherical  Hankel  function  of  the 

second  kind.  The  angle  function  normalization  N is  given  by 

un  w 


J^lS^h.n^dn 


2 I 

r*0,l 


• (r+2u) Id  (hlun) 

, rl(2r+2utl)  » 


where  r ■ 0,2,U,...  if  n-u  is  even  and  r»l,3,5,...  if  n-u  is  odd. 


Also, 


D(s,t,r,u,m) 


p£‘u(0)a(s,t,r,u,m)  + (-l)upj+u(0)a(8 


,t,r,-u,mj, 

(12) 


vhere 


a(a,t,r,u,m) 


^s+t-r 


f 

t (2s+l ) (2t+l ) (a-u)  I (t-a+u)  I (r+m) ! [ (a+t»r)/2  ] 1 
[ ( r+t-s ) /2 ] ! [ ( r+s-t ) /2 ] ! [ ( s+t-r ) / 2 ] ! ( s+t+r+1 ) I 


with 

(1*) 

(15) 

(16) 


(j)  - al/[bl  (a-b)  1 ] , 

wain  ■ 1/2  [Max  (r-s-t,s-r-t-2a*2u,t-s-r+2u)] , 
waax  ■ 1/2  [Min  (a+t-r,r+t-s-2a*-2u,r+s-t+2u)] , 


and 

p*(0)  . (-1  >(p-*)/2 

P ^(Kl),  ,£11, , • 


(17) 


since  p-q  = t-m±u  is  required  to  be  even.  The  quantity  is  equal 
to  1 if  u » 0 and  is  equal  to  2 if  u j 0. 

Expressing  p(r)  entirely  in  by  use  of  (8),  differentiating  with 
respect  to  in  order  to  apply  the  boundary  condition  on  given  by  (5), 
Multiplying  by  8^(h,n)cos  q+  and  integrating  over  n frost  -1  to  ♦!  and 
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over  4 from  0 to  2w  , the  ranges  of  orthogonality  for  these  functions 
one  obtains  the  following  set  of  equations: 


where  the  dot  indicates  the  derivative  with 


Similarly , application  of  the  boundary  condition  on  leads  to  the 
additional  set  of  equations: 


In  theory,  the  infinite  set  of  equations  given  by  (18)  plus  (20)  can  be 
•olved  for  the  unknowns  A>(  and  B^t  . In  practice,  the  summations  in  (18) 
and  (20)  are  truncated  at  1 ■ L,  and  the  remaining  T • (L+l)(L+2)  equations 
are  solved  for  A and  B , l ■ 0 to  L,  ■ ■ 0 to  t . The  choice  for  L 


depends  on  the  desired  accuracy  of  the  results,  the  normal  velocity 
distribution  vfi^),  and  the  acoustic  size  ka.  An  estimate  for  L is  given 
by  the  highest  order  that  must  be  retained  in  the  expansion  of  v(n^)  in 
*•*■»  ©f  spheroidal  angle  wave  functions.  The  truncated  equations  can  be 


rnmamtmm 


written  in  matrix  form 


* 

l M br  - g , q - 1 to  T , 


where  the  matrix  elements  depend  only  on  the  frequency  and  the 
geometry  of  the  two  disks , the  elements  b^  consist  of  the  unknown  coef- 
ficients and  , the  right  hand  side  elements  g^,  q ■ 1 to  T/2, 
depend  on  the  normal  velocity  distribution  v(n^) » and  the  elements 
g^»  q * T/2+1  to  T are  zero.  Equation  (21)  can  be  inverted  to  give 


br  "q-1  M^8q  ’ r " 1 t0  T ’ 


so  that  once  the  inverse  elements  have  been  calculated  for  a given 
frequency  and  geometry,  kml  and  can  be  obtained  for  any  velocity 
distribution  on  by  matrix  multiplication.  The  method  can  easily  be 
extended  to  include  a nonzero  velocity  distribution  on  by  modifying 

the  right  hand  side  of  (20)  to  give  nonzero  elements  g , q ■ T/2+1  to  T. 

4 

Consider  the  special  case  of  uniform  vibration  where  the  normal 
velocity  distribution  v(n^)  on  is  equal  to  a constant  over  the 
part  of  called  S+  , is  equal  to  a constant  -V^  over  the  part  of 
called  S , and  is  equal  to  zero  over  the  remainder  of  D^.  For  example, 
the  pulsating  disk  has  an  S+  region  covering  the  entire  surface  -lsn^sl  , 
while  the  oscillating  disk  has  both  an  S+  region  covering  Osn^sl  and  an 
8 region  covering  -lsn^,  <0  . It  is  still  assumed  that  is  stationary. 

Under  these  conditions  the  normalized  self  acoustic  radiation  Impedance 
ss  is  defined  by 
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quantities  rg  and  xg  are  the  nornalized  self  acoustic 
e and  reactance,  respectively.  Expressing  p(r)  entirely 
and  substituting  into  (23),  one  obtains 


depending  on  whether  l 


ouu.Lar.Ly,  the  force  produced  by  p(r)  on  that  p 
dirt  D2  vhicb  corresponds  to  the  vlbr.tlng  eree  on  D. 

end  used  to  ohteln  the  nomellsed  .utuel  eooustlc  rm, 
This  gives 
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0CV1A1  ^ P^n2’*2^dS^n2**2^ 


-ft  . 

I I A . [* 

i-0  Lja-0  “ n»0,l 


(27) 


ecVlAl 

+ BOlR01^h,0^Ioij  » 


CaUOnROn^h,O)l0n 


(28) 


where  the  integration  is 
S+  and  S". 


over  the  sane  range  in  r\  as  covered  by  n in 


When  both  disks  are  vibrating  uniformly  with  the  sane  velocity 
attribution  and  with  individual  velocity  snplitudes.  ^ «d  v2,  the 
total  nornalized  acoustic  radiation  impedances  ^ for  ^ and  ^ tor  D 

^ °btained  **  the  lnP®dances  z#  and  zffi  by  the  expression. 


*1  • *.  ♦ *»VV1  • 
*2  * *.  ♦ Vl/V2  • 


(29) 

(30) 


““  «,  and  .B  characterize  the  aconatic  radiation  lnp.d«c«  of  the  teo 
dirt  ayatem.  I„  the  renainder  of  the  paper  t,.  r.  , z,  and  ^n  . r*  , 
«ill  be  referred  to  aa  the  ..if  radiation  tnpedance.  realatance,  Jd 

reactance  and  the  mtuai  radiation  lapedance,  reai.taac.  w reactance, 
respectively. 
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NUMERICAL  RESULTS  AND  DISCUSSION 
a.  Introduction 


Calculation  of  the  self  and  mutual  radiation  impedances  by  use  of 
(25)  and  (28)  requires  numerical  values  for  the  oblate  spheroidal  radial 
wave  functions  R^^(h,0)  and  R^^(h,0)  and  their  first  derivatives  vith 
respect  to  £ and  for  the  integrals  1^^  and  . The  necessary  radial 
functions  were  calculated  using  OBRAD  [ 5 ] , a Fortran  computer  program 
written  in  double  precision  arithmetic  for  the  CDC  3800  computer  at  the 
U.  S.  Naval  Research  Laboratory*  OBRAD  is  one  of  several  computer  programs 
[8*7]  which  have  been  developed  recently  to  evaluate  spheroidal  wave 
functions  with  greater  accuracy  over  a wider  range  of  parameter  than 
previously  available.  Extensive  tables  [8,9]  of  both  oblate  and  prolate 
spheroidal  radial  wave  functions  and  their  first  derivatives  with  respect 
to  £ have  recently  been  published.  These  tables  contain  entries  for  values 
of  m » 0,1,2;  t ■ m(l)m+l»9  for  a wide  range  of  values  of  £ and  h.  The 
integrals  *0i  and  1^  were  evaluated  by  expanding  the  angle  functions  in 
a series  of  Legendre  functions  as  given  in  (If),  subsequently  expanding  the 
Legendre  functions  in  a series  in  cosrO  , where  6 « cos_1n  and  r is  an 
integer,  and  evaluating  in  closed  form  the  resulting  simple  integrals. 

b.  Pulsating  Disk  and  Stationary  Disk 


Consider  first  the  case  where  dish  is  pulsating  with  a normal 
velocity  V^,  and  disk  Dg  is  stationary.  Since  the  z ■ 0 plane  exclusive 
of  the  surface  of  D^  can  be  replaced  by  a stationary  surface  without 
affecting  the  acoustic  fields,  the  situation  is  equivalent  to  the  uniform 
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vibration  of  a circular  piston  in  a rigid  infinite  plane.  This  latter 
problem  has  been  solved  previously.  A detailed  review  by  Bouwkamp  is 
found  in  [10  ].  Extensive  tables  of  both  z and  z for  various  values  of 

8 B 

ka  and  kr^  calculated  from  expressions  given  by  Rayleigh  [ll]  and 
Pritchard  [12]  vere  published  recently  by  Berkovitz  and  Baler  [13].  The 
existence  of  these  independent  results  provides  an  accurate  check  on  the 
method  described  in  this  report.  Corresponding  values  for  z and  z vere 

8 B 

obtained  from  (25)  and  (28)  for  numerous  values  of  ka  and  kr^  end  compared 
with  the  results  of  Berkovitz  and  Baler.  The  agreement  vas  never  less 
than  3 significant  figures. 

c.  Oscillating  Disk  and  Stationary  Disk 

Consider  the  case  where  disk  is  oscillating  vith  a normal  velocity 
amplitude  V^,  and  disk  is  stationary.  The  self  radiation  impedance  of 
vas  calculated  from  (25).  It  vas  found  that  the  self  radiation  reactance 
xs  does  not  vary  significantly  from  that  of  a single  oscillating  disk  as 
the  spacing  kr12  is  decreased  from  over  200,  where  the  effect  of  the 
stationary  disk  is  negligible, to  2ka,  where  the  two  disks  are 

touching.  These  results  are  summarized  in  Table  1.  The  self  radiation 
resistance  rg,  however,  is  somewhat  sensitive  to  kr12.  Figure  2 shows  r^ 
as  a function  of  the  relative  spacing  r^/a  for  several  values  of  ka. 

Because  of  the  large  dynamic  range  of  the  values  for  r#,  some  of  the 
curves  have  been  scaled  as  indicated  in  the  figure.  Note  that  the  effect 
of  the  stationary  disk  is  significant  vhen  both  r^/a  and  ka  are  small, 
but  decreases  rapidly  as  either  becomes  large.  Vhen  r^/a  is  greater  than 
about  U.O  or  ka  is  greater  than  about  1.5*  the  effect  of  the  stationary 


104 


disk  is  negligible  so  that  is  very  nearly  equal  to  the  radiation 
resistance  of  a single  oscillating  disk.  Crane '[2]  investigated  the  self 
radiation  impedance  for  this  two  disk  problem.  His  calculations,  restricted 
to  ka  ■ 0.3  to  1.0  and  involving  approximations  of  uncertain  validity, 
also  shoved  that  the  stationary  disk  has  a negligible  effect  on  the  self 
radiation  reactance.  Hovever,  his  results  for  the  self  radiation  resist- 
ance, although  in  qualitative  agreement  with  Figure  2,  show  an  effect  of 
the  stationary  disk  that  is  many  times  larger  than  that  obtained  by  the 
present  method. 

The  mutual  impedance  between  the  oscillating  disk  and  the  stationary  disk  vas 

calculated  by  the  use  of  (28).  The  behavior  of  both  the  mutual  radiation 

resistance  r and  the  mutual  radiation  reactance  x as  a function  of  kr,  „ 
m m 12 

is  shown  in  Figures  3-6  for  several  values  of  ka.  It  is  seen  that  the 

effect  of  the  oscillating  disk  on  the  stationary  disk  tends  to  increase  as 

ka  is  increased  until  about  ka  ■ 2.5  when  it  then  begins  to  decrease. 

Mote  that  point  dipole  theory  predicts  that  r and  x will  be  identically 

n xb 

equal  to  zero  since  the  acoustic  pressure  vanishes  throughout  the  z plane. 

Thus  the  initial  increase  in  the  effect  as  ka  increases  from  small  values 
is  due  to  the  deviation  of  the  pressure  field  from  that  of  a point  dipole. 

Both  the  finite  size  of  the  oscillating  disk  and  scattering  from  the 
stationary  disk  are  contributing  factors.  As  ka  continues  to  increase, 
however,  the  competing  effects  of  increased  directivity  associated  with 
the  radiation  from  the  oscillating  disk  and  increased  phase  cancellation 
over  the  stationary  disk  eventually  dominate,  and  the  effect  begins  to 
decrease. 
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One-Sided  Vibrating  Disk  and  Stationary  Disk 


How  consider  the  case  where  disk  is  vibrating  uniformly  over  only 
its  top  half  with  a normal  velocity  amplitude  V^.  The  bottom  half  of 
and  all  of  D2  are  assumed  to  be  stationary.  Note  that  the  surface  normal 
velocity  distribution  for  this  problem  is  equal  to  one  half  of  the  sum  of 
the  surface  normal  velocity  distributions  for  the  previous  two  examples. 
Consequently,  the  acoustic  field  for  this  problem  can  be  obtained  as  one 
half  of  the  sum  of  the  acoustic  fields  for  the  separate  problems:  (a) 

oscillating  and  stationary  (b)  pulsating  and  stationary. 

This  example  of  linear  superposition,  knovn  as  the  Gutin  principle  [ lU ] 
results  from  the  linear  wave  equation  and  is  represented  in  the  linear 
nature  of  equations  such  as  (19)  and  (22).  As  mentioned  earlier,  the 
acoustic  radiation  from  (b)  is  equivalent  to  the  radiation  from  a uniformly 
vibrating  circular  piston  in  a rigid  infinite  plane  and  can  be  evaluated 
accurately  by  the  use  of  formulas  given  in  the  literature. 

The  self  radiation  impedance  of  was  calculated  for  several  values 
of  ka  and  kr^2  by  the  use  of  (25).  The  results  are  summarized  in  Table  2. 
Both  the  self  radiation  resistance  and  reactance  remained  virtually  unchanged 
from  the  corresponding  values  for  a single  disk  vibrating  from  one  side  as 
the  spacing  kr12  is  decreased  from  over  200  to  2ka.  Thus  the  effect  of  the 
stationary  disk  on  the  self  impedance  of  the  one-sided  vibrating  disk  is 
negligible  for  all  values  of  ka  and  kr^2«  Crane  [2]  presented  values  for 
r#  that  were  obtained  by  combining  his  values  for  the  self  radiation 
impedance  of  an  oscillating  disk  in  the  presence  of  a stationary  disk  with 
values  for  the  self  radiation  impedance  of  a circular  piston  in  a rigid 


106 


I 

I 

I 

! 

J 


L 


r 

i . 


r 

i. 


ii 


i 


i; 


E 


e 


E 


infinite  plane.  Note  that  the  radiation  impedance  of  the  circular  piston 
does  not  depend  on  kr12,  so  that  the  entire  variation  of  r#  vith  kr^2 
COB**  ^ro®  the  oscillating-dlsk  contribution.  Therefore  the  apparent 
inaccuracy  in  Crane's  values  for  the  oscillating-disk  contribution  results 
in  a variation  of  rg  vith  kr^  that  is  too  large  (up  to  20% ) . 

The  mutual  radiation  impedance  betveen  the  vibrating  half  of  and 
the  top  half  of  the  stationary  disk  was  calculated  by  the  use  of  (28). 
Values  of  r^  and  for  various  values  of  ka  are  plotted  versus  kr^2  in 
Figures  7 and  8 and  Figures  9 and  10,  respectively.  The  results  for  small 
**  *re  good  agreement  vith  calculations  using  monopole  theory  where 
is  replaced  by  a uniformly  pulsating  sphere  vith  an  equivalent  volume  accel- 
eration, and  Dg  is  considered  to  be  a single  point.  For  small,  increasing 
ka,  the  effect  of  the  one-sided  vibrating  disk  on  the  stationary  disk 
increases  as  the  volume  acceleration  of  also  increases  and  produces  a 
larger  acoustic  pressure  everywhere.  Hovever,  as  in  the  case  of  the 
o®ciHating  disk,  the  competing  factors  of  increased  directivity  and  phase 
cancellation  eventually  dominate  and  produce  a decrease  in  the  effect 

above  about  ka  ■ 2.0.  Hanish  [3]  calculated  values  of  r for  ka  * 0.1  to 

n 

1.0  by  use  of  the  linear  superposition  described  above.  He  used  Pritchard's 
formula  to  evaluate  accurately  the  contribution  from  the  circular  piston 
in  a rigid  infinite  plane  but  neglected  scattering  in  calculating  the 
contribution  from  the  oscillating-disk  component . Hovever,  since  the 
contribution  from  the  circular  piston  is  usually  far  greater  than  that 
from  the  oscillating-disk  component,  Hanish 's  results  are  usually  in  good 
agreement  vith  those  obtained  using  (28). 
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Self  radiation  reactance  of  an  oscillating  disk  in  the  presence  of  a 
stationary  disk 


ka 

kr^  - 2ka 

x. 

kr12  - 200 

x. 

0.1 

0.01*26 

0.0k25 

0.3 

0.130 

0.130 

0.5 

0.22k 

'0.223 

1.0 

0.51k 

0.515 

1.5 

0.888 

0.90k 

2.0 

1.035 

1.001 

2.5 

0.519 

0.533 

3.0 

0.225 

0.226 
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TABLE  2 


8*1 f radiation  impedance  of  a one-sided  vibrating  disk  in  the  presence 
of  a stationary  disk 
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FIGURE  CAPTIOHS 


Fig.  1 The  two  disk  geometry. 

Fig.  2 The  eelf  radiation  resistance  of  an  oscillating  disk  in  the 
presence  of  a stationary  disk  as  a function  of  the  relative 
spacing  r12/a. 

Fig.  3 The  mutual  radiation  resistance  between  an  oscillating  disk  and 
a stationary  disk  as  a function  of  the  separation  kr^. 

Fig.  k The  mutual  radiation  resistance  between  an  oscillating  disk  and 
a stationary  disk  as  a function  of  the  separation  kr^. 

Fig.  5 The  mutual  radiation  reactance  between  an  oscillating  disk  and  a 
stationary  disk  as  a function  of  the  separation  kr^g. 

Fig.  6 The  mutual  radiation  reactance  between  an  oscillating  disk  and  a 
stationary  disk  as  a function  of  the  separation  kr^  • 

Fig.  T The  mutual  radiation  resistance  between  a one-sided  vibrating 

disk  end  a stationary  disk  as  a function  of  the  separation  kr^  p . 

Fig.  8 The  mutual  radiation  resistance  between  a one-sided  vibrating 

disk  and  a stationary  disk  as  a function  of  the  separation  kr^. 

Fig.  9 The  mutual  radiation  reactance  between  a one-sided  vibrating  disk 
and  a stationary  disk  as  a function  of  the  separation  kr12. 

Fig.  10  The  mutual  radiation  reactance  between  a one-sided  vibrating  disk 
and  a stationary  disk  as  a function  of  the  separation  kr... 
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MODELING  OF  FINITE  ARRAYS  OF  COMPLIANT  TUBES 
R.P.  Radlinski 

Naval  Underwater  Systems  Center 

INTRODUCTION 

The  modeling  described  here  concerns  the  use  of  two  dimensional 
multibody  programs  to  predict  the  performance  of  a finite  width,  plane 
array  of  compliant  tubes.  The  current  plans  at  NUSC  are  to  imbed  at  least 
one  planar  array  of  plastic  compliant  tubes  in  a rubber  blanket  and  to 
test  this  composite  when  attached  to  a flexing  plate.  The  two-fold  task 
of  the  baffle  is  to  (1)  provide  a noise  free  environment  for  outboard 
mounted  hydrophones  by  reflecting  near  field  radiation  generated  by  flexural 
wave  back  toward  the  plate,  and  (2)  to  reduce  the  level  of  sound  radiated 
into  the  far  field.  In  this  portion  of  the  report,  we  will  be  investigating 
the  performance  of  an  isolated  planar  array  of  tubes. 

An  individual  compliant  tube  can  be  designed  to  have  its  first 
transverse  flexural  resonance  in  the  frequency  band  of  interest.  Thus, 
the  tubes  are  truly  compliant  only  in  a finite  bandwidth.  A schematic 
diagram  of  a single  compliant  tube  is  shown  in  Figure  1.  Mathematically, 
the  tube  is  often  taken  to  be  rectangular  or  elliptical  in  shape.  Toulis 
had  considerable  success  by  analyzing  the  performance  of  a tube  as  a bar 
clamped  at  both  ends.^ 

2 

The  first  tubes  designed  by  Toulis  were  metallic  and  used  in  an 
acoustic  lens.  Metal  tubes  were  also  used  as  reflectors  in  the  NUSC-North 
American  Autonetics  AFAR  system  now  in  operation  near  the  Azores  Islands. 
Plastic  tubes  have  been  used  in  a Johns-Hopkins  Applied  Physics  Laboratory- 
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SCHMATfC  DIAGRAM  OF  A SINGLE  COMPLIANT  TUBE 

FIGURE  1 
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Autonetics  collaboration  of  the  construction  of  a Luneburg  lens.  Compliant 
tube  technology  has  been  pursued  at  NUSC  largely  through  the  efforts  of 
Gerald  Brigham  who  was  formerly  associated  with  the  Autonetics  collaborations 
TWO  DIMENSIONAL  MODELING 


Let  us  now  consider  modeling  a finite  array  of  compliant  tubes.  At 
the  last  meeting  of  this  committee  at  NRL,  I described  the  use  of  a two 
dimensional  arbitrary  body  model  which  was  used  to  calculate  radiation 
patterns  and  impedances  for  line  arrays  near  finite  width,  soft  reflectors. 
A two  dimensional  analysis  was  possible  because  the  lengths  of  the  inter- 


i; 
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acting  bodies  were  long  compared  with  the  acoustic  wavelengths  of  interest. 
The  situation  is  the  same  with  compliant  tubes.  By  extending  the  program 
to  include  incident  plane  wave  radiation,  it  becomes  feasible  to  study 
reflection  or  transmission  from  a finite  array  of  tubes. 


Consider  an  array  of  compliant  tubes,  as  shown  in  Figure  2,  where 
the  individual  tube  is  considered  elliptical  in  shape. 


■ 


The  total  pressure  at  a point  t is  given  by  the  sum  of  the  incident  and 


where 


k)  is  the  source  strength 


k is  the  wavenumber 


zero  order  Hankel  function) 


In  this  formulation,  the  source  strength  integral  describes  the  scattered 
wave.  The  unknown  source  strength  is  determined  by  satisfying  the  prescribed 
boundary  conditions. 


The  tubes  will  be  assumed  to  be  resonance  devices  with  an 


mechanical  impedance  given  by 


where 


Frequency 


First  flexural  resonance  frequency  of  the  tube 
as  measured  in  air 


Surface  area  of  the  tube 


Two  dimensional  volume  of  the  tube 


Bulk  modulus  of  the  tube 


semi-minor  axis  of  the  e 


semi-major  axis  of  the  ellipse 


where  p is  the  density  of  the  medium  and 


3 

3n' 


is  the  normal  derivative  at  the  surface  of  a tube. 


A comparison  of  the  prediction  of  the  model  with  measurements  is 
shown  in  Figure  3 for  an  array  of  ten  steel  tubes.  The  non-dimensional 
semiaxes  of  these  tubes  were  kb  = 0.138tt  and  the  semiminor  axes  were 
kb  = 0.042  at  the  single  tube  resonance  frequency  as  measured  in  air.  The 
major  axis  of  a tube  was  aligned  parallel  to  the  direction  of  incident  plane 
wave  radiation  and  tlu.  center-to-center  spacing  was  kd  = 7 ii  at  the  above- 
frequency.  The  bulk  modulus  of  the  steel  tubes  were  calculated  and  measured 
to  be  about  1/11  that  of  water.  As  seen  from  the  measurement  in  the  plane 


of  the  tubes,  the  reactive  loading  effect  decreases  the  effective  resonance 


frequency  of  the  array  in  water  13  percent  with  respect  to  the  single  tube 
resonance  frequency  in  air.  The  predictions  from  the  arbitrary  body  program 
show  a frequency  shift  of  about  4 percent  as  compared  with  the  measurements, 
but  the  predicted  bandwidth  agrees  well  with  the  data. 

To  obtain  convergence  for  ten  tubes  with  the  above  program  requires  at 
least  five  minutes/frequency  of  computer  time  on  the  Univac  1108.  If  one 
wished  to  investigate  the  properties  of  larger  arrays  of  tubes  and  vary 
several  array  parameters,  such  an  analysis  becomes  quite  expensive.  A 
mathematically  simpler  approach  is  to  treat  the  tubes  as  acoustically  round 
but  retaining  the  same  mechanical  impedance  function  as  was  devised  for  the 
elliptical  shape.  Since  the  acoustic  wavelengths  of  interest  are  large 
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compared  with  the  dimensions  of  a tube,  this  approximation  should  not  be 
critical.  The  scattered  wave  can  then  be  described  by  series  of  cylindrical 
harmonic  wave  functions.  The  two  dimensional  volume  of  the  round  tubes  was 
taken  to  be  equal  that  of  the  corresponding  elliptical  tube. 

As  seen  from  Figure  4,  the  total  pressure  at  point  t becomes 


P * P + P 
t rinc  rscatt 


Pt=  e 


ikr  cos  0 


L °° 

* ? £ 

£*•1  q=-<x 


ZA  H (kr  )eiq0£ 

q q £ 


where  H (kr  ) is  the  qth  order  Hankel  function  and  A is  an  undetermined 
4 q 

series  coefficient. 


GEOMETRY  FOR  THE  ROUND  TUBE  ANALYSIS 

FIGURE  4 
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To  satisfy  the  boundary  conditions  at  the  Cth  cylinder,  rewrite 


W 


P.  in  terms  of  Bessel  functions  as 
inc 


P.  = e 
inc 


ikr 


cos0  . V'T  /,  \ iq(0„+4) 

°Jl 2^Jq(kr£)e  l 2' 


i. 


and  use  the  Addition  Theorem 


i 


iq0 . 00 

vkve  ■ T, 


J (kr  )H  (kr„,)em9^'(m  q)6*j 
m H m-q 


m= 


which  expresses  the  cylindrical  harmonics  with  respect  to  one  coordinate 
system  in  terms  of  a second  coordinate  system. 


The  predictions  from  the  round  tube  analysis  are  compared  in  Figure  5 
with  measurements  taken  in  the  center  of  the  ten  steel  tube  array.  The 
frequency  shift  of  8 percent  for  the  predicted  curve  with  respect  to  the 
measurements  is  more  pronounced  than  the  curve  obtained  with  the  elliptically 
moded  tubes  (Figure  3).  This  additional  frequenty  shift  is  in  part  due  to 
the  closer  effective  spacing  of  the  round  tubes  with  respect  to  the  equal 
area  squashed  elliptical  tubes.  As  reported  by  Toulis*”,  the  mutual  coupling 
in  a plane  area  of  compliant  tubes  causes  the  resonance  frequency  to  approach 
that  in  vacuum  as  the  spacing  between  the  tubes  decreases. 


CONCLUSION 

Recent  measurements  made  at  NUSC  with  closely  spaced  plastic  tubes 
will  be  used  to  determine  any  further  modifications  which  might  be  made  to 
the  above  models.  Comparisons  of  model  predictions  with  measurements  from 
an  area  of  ten  steel  tubes  indicate  that  (1)  the  impedance  function  used 
to  describe  the  tubes  is  reasonable  for  sparsely  spaced  arrays  and  (2) 
because  the  acoustic  wavelengths  are  large  compared  with  the  dimensions  of 
the  tube,  the  mathematical  shape  used  for  an  individual  tube  is  not  critical. 
The  models  described  here  are  only  the  first  step  toward  understanding  the 
array  performance  in  the  actual  situation  where  the  tubes  are  imbedded  in 
rubber  and  attached  to  a flexing  metal  surface. 
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